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ABSTRACT 



We present Chandra X-ray Observatory ACIS-S3 X-ray imaging observations and VUr/F0RS2 and Hubble Space Telescope optical 
observations of two low-density Galactic globular clusters; NGC 6366 and M55. We detect 16 X-ray sources with 0.5-6.0 keV lumi- 
nosities above Lx = 4 x 10'" erg s"' within the half-mass radius of M55, of which 8 or 9 are expected to be background sources, and 
5 within the half-mass radius of NGC 6366, of which 4 are expected to be background sources. Optical counterparts are identified 
for several X-ray sources in both clusters and from these we conclude that 3 of the X-ray sources in M55 and 2 or 3 of the X-ray 
sources in NGC 6366 are probably related to the cluster. Combining these results with those for other clusters, we find the best fit for 
a predicted number of X-ray sources in a globular cluster ji^ = 1.2T + l.lMh, where T is the collision number and Mh is (half of) 
the cluster mass, both normalized to the values for the globular cluster M4. Some sources tentatively classified as magnetically active 
binaries are more luminous in X-rays than the upper limit of Lx - O.OOlLboi of such binaries in the solar neighbourhood. Comparison 
with XMM and ROSAT observations lead us to conclude that the brightest X-ray source in M55, a dwarf nova, becomes fainter in 
X-rays during the optical outburst, in accordance with other dwarf novae. The brightest X-ray source in NGC 6366 is a point source 
surrounded by a slightly offset extended source. The absence of galaxies and Ho- emission in our optical observations argues against 
a cluster of galaxies and against a planetary nebula, and we suggest that the source may be an old nova. 

Key words. Globular clusters: individual (NGC 6366 and M55) 



1. Introduction 

All stars emit X-rays, but some emit more than others. Thanks 
to the Chandra X-ray Observatory the study of X-ray sources in 
globular clusters now includes sources down to luminosities of 
Lx ~ 10^'^"^°ergs"' in the 0.5-2.5 keV band. At these low lu- 
minosities, most sources are low-mass main-sequence stars that 
rotate rapidly, which in old stellar clusters is the case only when 
they have been spun up by tidal forces or by accretion in close 
binaries. At the high luminosity end, Lx ^ 10^^ ergs"', most 
sources are low-mass X-ray binaries in which a neutron star ac- 
cretes matter from a companion. At intermediate luminosities, 
most sources are cataclysmic variables (CVs). The clusters best 
studied at low X-ray luminosities are 47 Tuc (Heink e et al.|2005 
and references therein) and M4 (Bassa et al.^ 2004) . 

Whereas the binaries of main-sequence stars are primordial, 
i.e. formed as binaries when the component stars formed, the 
low-mass X-ray binaries in globular clusters are thought to be 
formed in close stellar encounters that bring a previously sin- 
gle neutron star into a binary (see the review by |Verbunt &] 
Lewin|[2006 and references therein). For cataclysmic variables 
both formation mechanisms, evolution from a piimordial binary 
and capture of a previously single white dwarf in a close en- 
counter, are viable, depending on the circumstances. The pro- 
genitor binary of a cataclysmic variable must be wide enough 



to allow the more massive star to evolve into a fairly big gi- 
ant, before it reaches its Roche lobe. Such wide binaries are de- 
stroyed ('ionized' ) in cores of g lobular clusters with densities 



po > 10"^ Mopc ^ (Davies 



1997b. 



One would naively expect that all cataclysmic variables in 
these cores were formed in stellar encounters. However, it has 
been pointed out that cataclysmic variables that evolved in the 
low-density outskirts of globular clusters can sink towards the 
core at late times. Thus, the population of cataclysmic variables 
in dense cores can be a mixture of locally produced products of 
stellar encounters and recently arrived products of the evolution 



of primordial binaries ( jlvanova et al. 2006 ^ see also Hurley et al. 
[20571 ). 

The number of stellar encounters in a globular cluster scales 
roughly with the collision number F = po'^r^^ /v, where po is the 
central density of the cluster, r^. the core radius, and v the central 
velocity dispersion. Through the virial theorem v o^r^^sfpQ and 



thus F oc 



„ 1.5„ 2 



(Verbunt & Hut 1987 Verbunt 2003). The 
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mass M of a cluster can be estimated from the total luminosity 
of a cluster, with use of a mass-to-light ratio appropriate for a 
cluster star population, and if necessary a bolometric correction. 
The number of X-ray sources in a globular cluster then may be 
expected to depend on the collision number (dynamical origin) 
and on the cluster mass (primordial origin); the form of the de- 
pendence depends on the type of binary, and therefore on the 
luminosity range that is studied. 



C. G. Bassa et al.: X-ray and optical observations of M55 and NGC 6366: evidence for primordial binaries 



Table 1. Parameters of M55 and NGC 6366 used in this pa- 
per; core radius r^, half-mass radius r^, distance d and red- 



dening Eb-v - Av/3.1 (taken from Harris 1996 version Feb 
2003). The interstellar column A^h is cornputed from A^h - 



\.19x\(fi^cm-^A 



V, see 



Predehl & Schmitt ( 1995 1. 



Cluster 


'•cO 


rhC) 


d (kpc) 


Eb-v 


A'h (cm-^) 


M55 
NGC 6366 


2.83 
1.83 


2.89 
2.63 


5.3 
3.6 


0.08 
0.71 


4.44 X 10-" 
3.94 X 10-' 



The brightest low-mass X-ray binaries, with Lx > 
lO^'^ergs"', are sufficiently rare that a cluster contains either 
none or one; in one case (M15) two. The probability that a cluster 
contains such a source scales with F ( jVerbunt & Hut[[T987| ). The 
fainter low-mass X-ray binaries, with 10^"^ ^ Lx $ 10^^ ergs ', 
are more numerous; the number of them in a cluster also scales 



with F peinke et al.j |2003) [Gendre et al] |2003, Pooley et al 



|2003| l. These results appear to confirm the theoretical idea that 
neutron stars enter binaries in globular clusters via stellar en- 
counters. Studies of globular clusters within other galaxies com- 
plicate the issue, because they indicate that the probability that a 
cluster contains a bright source depends on its metallicity, being 
higher for metal-rich clusters (e.g. Kundu et al.|2002] l, and fur- 
thermore scales with po'''', a shallower dependence on density 
than implied by a proportionality with F (Jordan et al. 2004 1. 

The numbers of less luminous sources, mostly cataclysmic 
variables - hitherto only observed in globular clusters of our 
Milky Way - appear to have a more shallow dependence on 
central density A^ oc po'' with y - 0.6 - 0.7 (Heinke et al. 
2003] [Pooley et al.] [2053] ). [Heinke et &\\ ((2006| find different 



values of F for bright and faint ca taclysmic variables , with the 



boundary near 10 erg s . Finally, Kong et al. ( 2006 1 show that 



the number of sources (2-5 only!) in the low-density globular 
cluster NGC 288, a mixture of cataclysmic variables and bina- 
ries of main-sequence stars, scales with the mass of this cluster 
Pooley & Hut ( 2006|l successfully describe the number A^ of X- 



ray sources in globular clusters with an equation of the form 
N ^aM + bT. 

There are many uncertainties in the theoretical description 
of the formation of X-ray sources in globular clusters. We men- 
tion only the retained fraction of neutron stars; and the details 
of the spiral-in process which changes the wide orbit of the pro- 
genitor binary into the close orbit of a cataclysmic variable. To 
help untangle these uncertainties empirically we have observed 
a variety of globular clusters both in X-rays and in the optical. 
We have selected these clusters to span a wide range in central 
density, core radius, and mass. In this paper we discuss our ob- 
servations of two clusters with relatively large core radii and low 
central densities: M55 and NGC 6366. Some parameters of both 
clusters are given in Tables[T]andl6] 

Both clusters were observed with the ROSAT PSPC, which 
detected one source in the core of each cluster Pohnston et al.| 
1996[l. The position of the source in M55 was more accurately 



determined with a ROSAT HRI observation |Verbunt[ ( [2001| l, 
and this source was later optically identified with a dwarf nova 
(Kaluzny et al. 2005 1. XMM detected five sources in the core 
of M55, one of them coincident with the ROSAT source ( |Webb[ 
[etaLi[2506l l. 

2. X-ray observations 

M55 was observed for 33.7 ks on 2004 May 1 1, and NGC 6366 
was observed for 22.0 ks on 2002 July 5, both with the Advanced 



CCD Imaging Spectrometer (ACIS) on the Chandra X-Ray 
Observatory with the telescope aim point on the back-side il- 
luminated S3 chip. The data were taken in timed-exposure mode 
with the standard integration time of 3.24 s per frame and teleme- 
tered to the ground in faint mode. 

Data reduction was performed usingthe CIAO 3.3 software 
provided by the Chandra X-ray Center] The data were repro- 
cessed using the CALDB 3.2.2 set of calibration files (gain 
maps, quantum efficiency, quantum efficiency uniformity, ef- 
fective area) including a new bad pixel list made with the 
acis_runJiotpix tool. The reprocessing was done without in- 
cluding the pixel randomization that is added during standard 
processing. This omission slightly improves the point spread 
function. The data were filtered using the standard ASCA grades 
(0, 2, 3, 4, and 6) and excluding both bad pixels and software- 
flagged cosmic ray events. Intervals of strong background flaring 
were searched for, but none were found. The extraction of counts 
and spectra and the generat ion of response fi les was accom- 
plished with ACIS Extracfl proos et al. [[20021 version 3.107), 
which calls many standard CIAO routines. 

2. 1 . Source detection 

The CIAO wavelet-based wavdetect tool was employed for 
source detection in both the 0.5-6.0 and 0.3-10.0 keV bands. 
We detected 29 sources on the entire S3 chip in the M55 ob- 
servation and 12 in NGC 6366 observation. Of these sources, 
15 (14) lie within the 2'89 half-mass radius (2'83 core radius) 
of M55, and 3 (3) lie within the 2!63 half-mass radius (1'83 
core radius) of NGC 6366. We also searched part of the adja- 
cent S4 CCD in the M55 observation since part of the half-mass 
region fell on this chip, but no sources were detected in this area. 
Furthermore, we examined adaptively-smoothed images (made 
from the CIAO csmooth tool) for significant point sources. In 
the M55 observation, we found two additional possible sources 
within the core radius. In the NGC 6366 observation, we also 
found two additional possible sources within the half-mass ra- 
dius, one of which was within the core radius. Fig. [Tlshows the 
location of the detected sources with respect to the cluster center 
and the core and half-mass radius. 

All sources are consistent with being point sources, except 
for CXI in NGC 6366, which appears to extend to a radius of 
about 20" (Fig.|2|. Visual inspection suggests the superposition 
of a point source and extended emission. To assess their spatial 
coincidence, we fit two gaussians plus a constant background 
to an image of the region. One gaussian approximates the point 
spread function, and the other is broader and describes the ex- 
tended emission. Our best fit model gives source counts in rough 



agreement with the methods described in § 2.2 and given in 
Table[3| according to the model, the point source has 17.4 net 
counts, and the extended source has 186.1 net counts. The two 
gaussians are offset from each other by 2'.'4 + I'.'l. 

We will give two estimates of how many sources are related 
to each cluster For the first estimate, we assume that all sources 
outside the half-mass radius are fore- or background sources, 
A^o = 14 for M55 and 9 for NGC 6366. Taking the size of the 
S3 detector as 8!4 x 8f4 = 70.6 sq.arcmin, the expected number 
of fore- or background sources within the half-mass radius rj,, 
expressed in arcmin, follows as A^e = N^nr]^ I{1Q.6 - nr^^). This 
gives A^e = 8.3 for M55 and 4.0 for NGC 6366. 

' http://asc.harvard.edu 

- http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 
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Fig.l. 10' X 10' Digitized Sky Survey R-hand images of M55 (left) and NGC 6366 (right). The positions of the Chandra X-ray 
sources are depicted with circles and numbered according to their CX designation (see Tablel2]and[3 i. The two large rectangles on 
each image denotes the field-of-view of the two F0RS2 chips, while the smaller slanted rectangles denote the field-of-view of the 
Hubble Space Telescope ACS/WFC chips. 



The second estimate may be found from the log A^-log S re- 
lationships of Giacconi et al. (2001 1; at our detection limits, this 
leads to an expected number of background sources within the 
half-mass radius of 6-9 for M55 and 3-4 for NGC 6366. These 
numbers are wholly compatible with our first estimate. In M55 
we detect 17 sources within the core, of which we estimate 8-11 
cluster members. The probability to find 14 or more sources for 
expected 8.3 is about 4%; so it is likely that at least 4 sources 
are members. In the case of NGC 6366 we estimate that 1 of the 
4 sources within the half-mass radius is a cluster member, but 
there is a sizable probability (37%) that all sources are fore- or 
background. 

We repeat our first estimate for sources with Lx > 4 x 



10^" ergs ' [0.5-6.0 keV], the limit used in the ,Pooley et al. 



( 2003| l study. Above the corresponding flux Umit of 1.2 x 
0"" ergs"' cm"^ for M55, we detected 14 sources outside the 
half-mass radius of M55, and expect 8.3 sources unrelated to the 
cluster within the half-mass radius. 16 sources within the half- 
mass radius are observed above the flux Hmit. Above the flux 
limit for NGC 6366 of 2.6 X 10"'^ ergs"' cm'^, we observe 9 
sources outside the half-mass radius, which predicts 4.0 unre- 
lated sources within the core, compared to the observed number 
of 5. The probability to observe 14 or more sources when 8.3 are 
expected is about 4%. Thus from the source numbers alone we 
have marginal evidence for sources above the luminosity limit 
used by Pooley et al.| (|2003") related to M55; and no evidence for 
such sources related to NGC 6366. 

In Tablesl2]andl3]we list the properties of the X-ray sources 
detected by Chandra in M55 and NGC 6366, respectively. 

2.2. Count Rates 

We extracted source counts in the following bands: 0.5-1.5 keV 
(Xsoft), 0.5^.5 keV (X^ed), and 1. 5-6.0 keV (Xh^-d). The de- 
tected count rate was corrected for background, exposure vari- 



ations, and foreground photoelectric absorption. We make these 
corrections in order to produce a X-ray color-magnitude diagram 
(CMD) that can be compared to the X-ray CMDs that have re- 
sulted from Chandra observations of other globular clusters. In 
addition, however, attention must be paid to differences in detec- 
tor responses and, of course, exposure times and distances. 

The background count rates in each band were estimated 
from source-free regions on the S3 chip outside the half- 
mass radii. The density of background counts in each band 
for the M55 observation is 0.0050 counts pixel"' in Xjott, 
0.0103 counts pixel"' in Xmed, and 0.0072 counts pixel"' in Xhard- 
For the NGC 6366 observation, the background densities are 
0.0030 counts pixel"' in Xsos, 0.0063 counts pixel"' in Xmed, and 
0.0056 counts pixel"' in Xhard- The background count rates in the 
cores may be somewhat higher, but even factors of a few greater 
than this estimate have negligible effects on our analysis. 

The exposure variations among sources were at the 15% 
level or less in both observations. To account for these variations 
in exposure, we applied multiplicative corrections based on the 
ratio of the average effective area of the detector at the location of 
a source in each of the three bands to that in the same band of the 
source which had the highest average exposure (CXll in M55 
and CX2 in NGC 6366). The individual effective area curves for 
the sources were made using the CIAO tool mkarf. The aver- 
age effective area of the detector at the location of CXI 1 in M55 
in each of the bands was 494 cm^ (^soft), 449 cm^ (^med), and 
381 cm^ (Xhard); for CX2 in NGC 6366, the areas were 509 cm^ 
(X,oft), 447 cm2 (X^ed), and 373 cm^ (Xhard). 

While the previous corrections were relatively minor (at the 
few percent level or less), the correction for photoelectric ab- 
sorption is appreciable for NGC 6366 (less so for M55). We 
investigated the effects of absorption by the column densities 
given in Tablefl] on three characteristic spectra: a 3 keV ther- 
mal bremsstrahlung, a 0.3 keV blackbody plus power law with 
photon index of 2, and a power law with a photon index of 2. 
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Table 2. Chandra X-ray sources detected in our observation of M55. The positions of the X-ray sources have been coiTected for the 
bore-sight correction of -0'.'270 in right ascension and +0'.'080 in declination. Positional uncertainties are given in parentheses and 
refer to the last quoted digit and are the centroiding uncertainties given by ACIS Extract. They do not include the uncertainties in 
the bore-sight correction (Qf.'QAA in right ascension and Of.' 045 in declination). The X-ray bands are 0.5-1.5 keV (Xjoft), 0.5^.5 keV 
(^med) and 1.5-6.0 keV (Xhard)- The first 17 sources are located within the half-mass radius of this cluster and are ordered on 0.5- 
6.0 keV countrate (Xjoft +^hard)- The remaining sources are located outside the half-mass radius and are ordered on right ascension. 
Some of the X-ray sources are also detected by ROSAT ( John ston etlE] [19961 |Verbunt[ [200 1 p and XMM ( [Webb et al.[[2006| . These 



sources are denoted by R and X. CVl is cataclysmic variable found by 



Kaluznyetal. (20051. 



ID 


R.A. 


Decl. 


Counts 


(Detected/Cc 


rrected) 


/x (0.5-2.5 keV) 


/x (2.5-6.0 keV) 




(J2000) 


(J2000) 


^soft 


'^med 


^haid 


(erg s"' cm"^) 


(ergs-'cm-^) 


CXI (R9/X30/CV1) 


19H0"08^;593(3) 


-30°58'52'.'08(4) 


60/74.3 


102/123.5 


45/47.9 


1.2 X 10-'^ 


1.4 X 10-'" 


CX2 (X42) 


19''39"'57!877(3) 


-30°56'11'.'67(4) 


27/32.8 


44/51.4 


18/18.5 


5.0 X 10-'5 


5.8 X 10-" 


CX3 (XI 2) 


19h40"'04»720(3) 


-30°56'47'.'09(4) 


14/17.0 


33/38.6 


22/22.8 


3.4 X 10-'5 


7.9 X 10-" 


CX4 (X13) 


19H0"'09il92(5) 


-30°56'58'.'54(7) 


4/4.9 


21/24.7 


19/19.8 


1.9 X 10-'5 


7.2 X 10-" 


CX5 


19H0"08^;846(6) 


-30°56'20'.'34(8) 


10/11.8 


15/17.5 


7/7.5 


1.9 X 10-'5 


2.1 X 10-" 


CX6 


19''39"'53^;i08(ll) 


-31°00'14'.'99(15) 


5/6.0 


12/13.8 


8/8.2 


1.3 X 10-'5 


2.3 X 10-" 


CX7 


19''39'"5P.190(9) 


-30°59'01'.'01(11) 


10/11.8 


13/14.7 


3/3.1 


1.7 X 10-'5 


1.1 X 10-" 


CX8 


19''40"'03'>;279(8) 


-30°59'05'.'61(11) 


5/6.2 


10/12.1 


6/6.5 


1.2 X 10-'5 


1.6 X 10-" 


CX9 


19*'39'"54^;920(7) 


-30°57'18'.'01(8) 


0/0.0 


8/9.3 


10/10.2 


6.1 X 10-'« 


3.6 X 10-" 


CXIO 


19h4()m()4s932(ii) 


-30°59'37'.'68(15) 


3/3.7 


9/10.9 


6/6.5 


9.2 X 10-'" 


1.7 X 10-" 


CXll 


19H0"'0a;235(8) 


-30°57'39'.'85(10) 


3/3.3 


7/7.3 


4/3.6 


8.1 X 10-'« 


6.4 X 10-'" 


CX12 


19H0"'0H887(8) 


-30°56'19'.'48(10) 


2/2.4 


5/5.9 


4/4.2 


5.8 X 10-'" 


1.2 X 10-" 


CX13 


19*'39"58^;668(12) 


-30°59'22'.'01(16) 


0/0.0 


5/6.3 


6/7.0 


4.8 X 10-'« 


2.4 X 10-" 


CX14 


19''40'''0a;479(14) 


-30°59'38'.'08(18) 


0/0.0 


5/5.4 


5/5.1 


4.3 X 10-'« 


1.4 X 10-'^ 


CX15 


19''39"59!793(10) 


-30°55'26'.' 12(13) 


4/5.4 


5/6.3 


1/0.8 


6.4 X 10-"5 


5.8 X 10-'" 


CX16 


19H0"'06!829(10) 


-30°57'12'.'29(13) 


5/6.4 


5/5.9 


0/0.0 


7.6 X 10-'" 


2.0 X 10-" 


CX17 


19''40"'08'!214(12) 


-30°57'15'.'47(15) 


3/3.7 


4/4.6 


1/0.8 


5.7 X 10-'" 


6.9 X 10-'" 


CX18(R6/X9) 


19''39"'42^;i30(4) 


-30°55'45'.'48(5) 


71/101.4 


137/170.8 


72/75.8 


1.8 X 10-'* 


2.0 X 10-'* 


CX19 


19h39m44sg42(i3) 


-30°58'47'.'54(17) 


5/6.4 


11/12.8 


6/6.0 


1.2 X 10-'5 


2.0 X 10-" 


CX20 (XI 4) 


19^39"45^;876(4) 


-30°57'4a'03(4) 


29/38.5 


92/112.6 


78/82.9 


9.7 X 10-'5 


2.9 X 10-'* 


CX21 


19''39"'45^;989(19) 


-31°01'31'.'65(25) 


4/5.2 


14/15.2 


13/12.0 


1.5 X 10-'5 


4.7 X 10-" 


CX22 (XI 9) 


19*'39'M6!573(12) 


-3r02'26'.'47(15) 


33/44.2 


65/77.3 


34/33.9 


8.1 X 10-'5 


1.2 X 10-'* 


CX23 


19''39'M7!994(11) 


-30°55'24'.'12(14) 


8/23.8 


11/28.6 


3/6.5 


3.5 X 10-'5 


1.2 X 10-'5 


CX24 (R13/X17) 


19^39"55^;778(4) 


-3r02'04'.'98(5) 


207/264.9 


304/368.1 


106/112.5 


3.9 X 10-'* 


2.6 X 10-'* 


CX25 (X21) 


19''39"'57!005(15) 


-31°03'07'.'82(19) 


45/56.2 


48/56.0 


3/2.0 


7.1 X 10-'5 


4.8 X 10-'" 


CX26 


19H0"'0H219(18) 


-3r01'02'.'60(24) 


0/0.0 


3/3.0 


7/7.1 


6.9 X 10-'" 


3.9 X 10-" 


CX27 


19H0"'03^;284(24) 


-31°01'52'.'73(31) 


0/0.0 


8/9.4 


9/8.8 


6.4 X 10-'" 


3.1 X 10-" 


CX28 (X20) 


19H0"'03^;475(14) 


-3r02'42'.'40(18) 


25/30.0 


37/42.7 


13/12.6 


4.7 X 10-'5 


4.2 X 10-" 


CX29 (X45) 


19''40"'13'i402(15) 


-30°59'45'.'22(19) 


9/11.0 


11/13.1 


2/2.2 


1.6 X 10-'5 


7.7 X 10-'" 


CX30 


19H0"n6!842(13) 


-30°57'41'.'89(17) 


4/5.5 


10/12.3 


6/6.0 


1.0 X 10-'5 


3.0 X 10-" 


CX31 


19''40"U 8^:981(20) 


-31°00'13'.'35(26) 


2/2.2 


12/14.5 


12/12.1 


1.3 X 10-'"* 


5.0 X 10-" 



The effects were most prominent in the Xsott band. Averaging 
the results of each spectrum in each band, we use the following 
correction factors for NGC 6366: 3.72 (X^ott), 2.47 (X^ed), and 
1-17 (Xhard). For M55, the factors are 1.21 (Xsoft), 1-16 (X^ed), 
and 1.02 (Xhaid)- Table[2]and 3] list both the observed and fully 
corrected counts in each band. The effect of the absorption cor- 
rection on the X-ray CMD (Fig. [3} is a uniform shift of the 
NGC 6366 sources by 0.39 units on the left axis and 0.50 units 
on the bottom axis and a uniform shift of the M55 sources by 
0.06 units on the left axis and 0.07 units on the bottom axis. The 
bottom and left axes give the X-ray color and magnitude with- 
out this shift (they do, however, include the small corrections for 
background subtraction and exposure variations). 



From the best fit spectra, we calculated the unabsorbed fluxes, 
given in Tables[2][3] 

For the point source NGC 6366 CXI a, the fluxes quoted in 
Table[3]are from a fit of an absorbed power law to the unbinned 
spectrum, for which we use Cash statistics, and fix the absorp- 
tion to that of NGC 6366. This gave a power-law with photon 
index 3.5 ±0.6. For the extended source NGC 6366 CXlb we ex- 
tracted the spectrum from within a 40" radius centered on CXI a. 
We subtracted the spectrum of the point source CXI a, and used 
a nearby, source-free, region with 1' radius to estimate the back- 
ground. This gave a bremsstrahlung temperature of the extended 
component of kT 
in Table[3] 



2.0^^7 keV, and unabsorbed fluxes as listed 



2.3. Spectral Fitting 

We fit all sources with absorbed power-law spectral models in 
Sherpa ( [Freeman et al. 2001) using |Cash[ ( [ 197 9) statistics. We 
fixed the column density to the value given in Table[T] with only 
the power law photon index and normalization allowed to vary. 



3. Optical observations 

Optical observations of M55 and NGC 6366 were obtained with 
the FORS2 instrument at the Unit Telescope 1 (UTl) of the ESQ 
VLT in April and May 2005. Both globular clusters were ob- 
served in three filters, B, R and Ha, with exposure times chosen 
to maximize the dynamic range. Table[4[ provides a condensed 
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Table 3. Chandra X-ray sources detected in our observation of NGC 6366. The celestial positions of the X-ray sources have been 
corrected for the bore-sight correction of +(y.'Q3Q in right ascension and -0'.'256 in declination. Positional uncertainties are given 
in parenthesis and refer to the last quoted digit and are the centroiding uncertainties given by ACIS Extract. They do not include 
the uncertainties in the bore-sight correction ((X'OQl in right ascension and 0'.'082 in declination). The X-ray bands are as defined 
in Table|2l The first 5 sources are located within the half-mass radius of this cluster and are ordered on 0.5-6.0 keV countrate 
(Xsoft + ^hard)- The remaining sources are located outside the half-mass radius and are ordered on right ascension. CXI was also 
detected by ROSAT Johnston et aLl[T996HVerbunt|[200T] l. 



ID 


R.A. 


Decl. 


Counts (Detected/Corrected) 


/x (0.5-2.5 keV) 


fx (2.5-6.0 keV) 




(J2000) 


(J2000) 


^soft 


^racd 


^hard 


(erg s ' cm 2) 


(ergs-' cm-^) 


CXla 


17''27"4Z892(7) 


-05°05'05'.'28(11) 


10/33.2 


15/32.1 


5/4.8 


8.5 X 10-'' 


7.0 X 10-"- 


CXlb (R4) 


17i>27"'4Z78(7) 


-05°05'04'.'7(1.1) 


87/190.0 


185/280.9 


111/79.2 


4.6 X 10-'" 


4.9 X 10-'-* 


CX2 


17''27"'3S682(7) 


-05°04'47'.'63(11) 


2/7.5 


6/15.0 


5/5.9 


1.4 X 10-'^ 


2.8 X 10-'^ 


CX3 


17''27"48H70(19) 


-05°02'58'.'69(28) 


2/7.7 


5/11.5 


3/2.8 


1.3 X 10-'^ 


2.4 X 10-'^ 


CX4 


17''27"'4Z044(12) 


-05°03'55'.' 12(18) 


0/0.0 


4/10.8 


4/5.2 


5.5 X 10-'* 


3.1 X 10'^ 


CX5 


17h27"44U37(12) 


-05°04'27'.'20(17) 


3/12.1 


5/13.4 


2/2.6 


2.7 X 10-'^ 


3.6 X 10-"5 


CX6 


17''27"27:384(8) 


-05°04'4a'42(12) 


5/20.1 


15/38.6 


10/12.0 


4.6 X 10-'5 


3.7 X 10-'^ 


CX7 


17''27"'2S726(16) 


-05°03'0a'34(24) 


2/7.7 


6/15.0 


6/7.2 


1.7 X 10-'5 


3.6 X 10-'^ 


CX8 


17''27"29:883(6) 


-05°05'24'.'56(8) 


15/54.7 


17/40.9 


2/2.1 


1.5 X lO-'-t 


1.7 X 10-'« 


CX9 


171,270,32307(7) 


-05°06'56'.'56(10) 


2/6.7 


8/19.9 


6/7.1 


2.2 X 10-'5 


2.2 X 10-'5 


CXIO 


17''27"'36^:643(17) 


-05° 02' 31'.' 24(26) 


1/4.0 


2/5.4 


2/2.6 


1.8 X 10-'5 


5.1 X 10-'5 


CXll 


17''27"'4a:i38(7) 


-05°01'13'.'76(10) 


36/141.8 


80/207.4 


45/55.2 


2.9 X 10-'" 


1.6x10'" 


CX12 


17''27"4ff:533(9) 


-05°08'35'.'70(14) 


0/0.0 


5/12.7 


6/6.9 


7.0 X 10-"^ 


5.2 X 10-'^ 


CX13 


17''27"'5ff:335(14) 


-05°07'ia'27(21) 


4/15.1 


6/14.7 


2/2.1 


5.3 X 10-'5 


8.3 X lO-''' 


CX14 


17''27"'5Z830(11) 


-05°02'56'.'05(16) 


16/62.0 


29/71.4 


13/14.6 


1.0 X lO-"* 


5.5 X 10-'5 
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Fig. 2. (left panel) The adaptively smoothed brightness distribution of NGC 6366 CXI shows a point source and an extended source. 
The solid contour contains 90% of the flux in the Chandra point spread function (at the detector location of the source); the dashed 
contours indicate the extraction region for the counts of the extended source, (middle panel) Contours from the adaptively smoothed 
brightness distribution of CXI overlaid on the 15 minute Ha images. The white diagonal strip is due to the gap between the FORS2 
chips, see Fig. [T] There is no evidence for any extended Ha emission, (right panel) The cumulative energy distributions of the point 
source (large-stepped distribution) and the extended source. A Kolmogorov Smirnov test gives a probability of 1% that both energy 
distributions are the same, and we conclude that the extended source is harder. 



log of the observations. F0RS2 is a mosaic of two 2kx4k chips 
with a pixel scale of (7/126 pix-' . For the majority of the obser- 
vations, we used 2x2 on-chip binning, providing a pixel scale 
of (X' 252 pix-', except when the seeing was below 0'.'6, when no 
binning was applied. 

The images were reduced using the Munich Image Data 
Analysis System (MIDAS). All images were bias-subtracted and 
flatfielded using twilight flats. Next, we grouped the images shar- 
ing the same combination of chip, filter and exposure time. The 
images in each group were aligned using integer pixel offsets 
and co-added to remove cosmic rays and increase the signal-to- 
noise. Hence, for each cluster we obtained 18 separate stacked 
images, 9 for each chip. 



Both M55 and NGC 6366 are part of a survey of globular 
clusters with ACS/WFC on board the Hubble Space Telescope. 
M55 was observed for 284 s in a V-band filter (F6(86W) and 384 s 
in an /-band filter (F814W). For NGC 6366 the exposure times 
with 570 s in both F606W and F814W. Compared to FORS2, 
ACS/WFC has a smaller field-of-view (about 3!4 x 3!4) and con- 
tains only a few of the X-ray sources. For M55, CX2, CX3, CX8, 
CX9, CXll, CX12, CX13 and CX16 are coincident with the 
202" X 202" field of view, while for NGC 6366 CXI through 
CX5 are coincident. 
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Table 4. A log of the optical observations. Here, cr denotes the 
seeing and secz the airmass. 



Date & Time (UT) 


Exp. (S, R, Ha) 


cr(") 


secz 


NGC 6366 








April 8 07 


40-09 


20 


4x (5 m, 2 m, 15 m) 


0.6-0.7 


1.06-1.13 


April 9 07 


22-08 


10 


2x(5m, 2 m, 15 m) 


0.7-1.0 


1.09-1.16 


08 


13-08 


24 


6x(3s, Is, 9 s) 


0.7-1.1 


1.07-1.08 


08 


25-08 


47 


6x(30s, 12s, 90s) 


0.6-0.8 


1.06-1.07 


M55 








April 10 08:01-09:51 


4x(5m, 2 m, 15 m)" 


0.4-0.6 


1.06-1.33 


May 6 06:37-06:49 


5x(3s, Is, 9 s) 


0.6-0.8 


1.22-1.26 


06:50-07:05 


4x(30s, 12s, 90s) 


0.7-0.9 


1.18-1.22 


May 12 07:07-08:00 


2x(5m, 2 m, 15 m)'' 


0.4-0.5 


1.04-1.11 



i°g x'ioft/x-h 



' These observations were obtained with 1x1 binning. 



I°9 X soft/X hard 
0.5 



°X^. 





1 1 i 1 1 i 1 1 i 




1 


1 


■ 




NGC 6366 




lb 




- 








K^- 




■ 










i,^° 


- 




7 


e-^ 


„2 




o' 


: 




O^ 




, 1 . , , . 1 , . 


1 


, . 1 





. 



-0.5 



0.5 



l°9 x° ,,/x"; 



Fig. 3. X-ray color-magnitude diagram for the X-ray sources de- 
tected in the half-mass radius of M55 (top panel) and NGC 6366 
(bottom panel). The X-ray color is defined as the logarithm of the 
ratio of Xsoft (0.5-1.5keV) counts to Xh-^-d (1.5-6.0 keV) counts. 
The X-ray magnitude is the logarithm of X,ned (0.5^.5 keV) 
counts. The bottom and left axes provide the observed absorbed 
color and magnitude scales (^). 



3.1. Photometry 

We used DAOPHOT II, running inside MIDAS, for the pho- 
tometry of the stacked FORS2 images. We followed the recom- 
mendations of Stetson ( 1987') and obtained instrumental mag- 
nitudes using point-spread-function (PSF) fitting. The PSF for 
each stacked image was determined in an iterative manner First, 
we used an analytic PSF to remove stars in the proximity of those 
stars used to determine the PSF. This subsequently improved the 
PSF for the next iteration. We found that allowing the final PSF 
to quadratically depend on the position on the chip decreased the 
fitting residuals. We used the final PSF to determine instrumental 
magnitudes of all stars above a 5cr detection threshold. The re- 
sulting star-subtracted image was then searched for objects that 
were missed in the first iteration. These stars were added to the 
star list after which the PSF-fitting process was repeated. 

In order to match the stars found on the separate stacked 
images, we employed the following method. First, we deter- 
mined approximate coordinate transformations between the sep- 
arate images taken with the same chip of each cluster to create 
a single list containing all stars. Next, we removed multiple en- 
tries of the same star from the list by removing stars located 
with 3 pixels of another star The resulting master star list was 



then matched against the star lists for each of the separate im- 
ages, again using the approximate coordinate transformations. 
We determined improved, 6 parameter, coordinate transforma- 
tions from the matched master list and repeated the process of 
creating a master list and matching stars. 

Finally, the instrumental magnitudes for the separate chips 
were calibrated against photometric standards in the Mark A (for 
M55) and the Rubin 152 fields (for NGC 6366), using calibrated 
magnitudes from Stetson (2000) and fitting for zero-point and 
colour coefficients. Extinction coefficients in B and R were taken 
from the FORS2 webpage. The final rms residuals of the calibra- 
tions were 0.03 mag in B and 0.04 mag in R for the calibration of 
the NGC 6366 observations and 0.06 mag in B and 0.04 mag in 
R for the M55 observations. The Ha magnitudes were calibrated 
by adopting R-Ha - for the bulk of the stars. Fig. HI shows 
colour-magnitude diagrams (CMDs) for both globular clusters. 

For both clusters, we estimate 5cr limiting magnitudes of B = 
26.2 and Ho- - 24.8. The /?-band magnitudes are limited through 
the calibration, which relies on the instrumental B - R colors 
and which is not available if a star is not detected in the B-band 
images. Hence, the 7?-band limit depends on the B - R color 
through R = 26.2 - (B - R). 

For the ACS/WFC observations we used the photometry pre- 
sented in Anderson et al.^ (2008 ), who use highly specialized 
methods to determine accurate photometry of ACS/WFC im- 
ages. For our purposes, we have used the photometry that was 
transferred from the ACS/WFC filters into ground-based V and 
/-band magnitudes. Fig.lS] shows the CMDs for both globular 
clusters. 



3.2. Astrometry 

We aim to tie the stacked /?-band images for each globular 
cluster and each CCD to the International Celestial Reference 
System (ICRS), to astrometrically calibrate our photometric cat- 
alog and to allow for a comparison between the X-ray and optical 
frames. For each globular cluster and each CCD, we calibrated 
the stacked 1 s iJ-band images against UCAC2 astrometric stan- 
dards (Zacharias et al. 2004| l. Even though the stellar density in 
the cores of these globular clusters is higher than a typical stel- 
lar field, our astrometry did not suffer from severe blending of 
stars, and for each image only 3 to 6 stars of the 40 to 80 that 
overlapped with each image were rejected. Fitting for zero-point 
position, scale and position angle, we obtained astrometric solu- 
tions for each image that had root-mean-square (rms) residuals 
of (y.'04 to (y.'Ql in right ascension and declination. These solu- 
tions were transferred to the deeper /?-band images using the cal- 
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Fig. 4. Optical color-magnitude diagrams for M55 and NGC 6366 constructed from the FORS2 data. Candidate optical counterparts 
to the Chandra X-ray sources are indicated and numbered. 



ibrated positions of stars on the 1 s images and used to calibrate 
the photometric catalog for each globular cluster Typically, each 
solution used some 1300 to 1900 stars, resulting in rms residuals 
ofa.'Ol toa.'03. 

The drizzled ACS/WFC F606W images were calibrated 
against the F0RS2 observations. Due to the large number of 
stars (500 to 700) coincident between the FORS2 and ACS/WFC 
images, the astrometric calibration has rms uncertainties of less 
than 0'.'02 in right ascension and declination. 

To place the Chandra frame onto the ICRS we use identifi- 
cations based on the Chandra X-ray sources alone. We identify 
M55 CXI with the dwarf nova CVl discovered by Kaluzny et al. 



(200 5 ). From the position of the optical counterpart, using the 
astrometry presented above, we determine a Chandra bore-sight 
correction of Aa = -a.'270 ± a'044 and A6 = -ntt'OSO + 0'.'045. 
For NGC 6366, we find that CXO J172729.9-050524 coincides 
with the bright star (V = 10.7) BD-04°4280 (PPM 706759; 
|Roser & Bastian||1988[ ). This star is saturated in our observa- 
tions, but an accurate position is provided in the 2nd version of 
the U SNO CCD Astrograph Catalog (UCAC2; [Zacharias et"aL 
2004). The bore-sight correction is Aa = +a'030 + 0'.'091 and 



Ad — -0'.'256 ± 0'.'082. These offsets put the X-ray positions 
onto the International Celestial Reference System (ICRS). We 
note that the uncertainty in both corrections is dominated by 
the uncertainty in the X-ray position. Both corrections, however, 
are within the 0'.'6 (90% confidence) accuracy in the pointing of 
Chandra ( |Aldcroft et aL||2000| ). 

With the optical and X-ray astrometry, the final uncertainty 
on the position of a X-ray source is the quadratic sum of the 
centroiding uncertainty of the X-ray source (those tabulated in 
Tablel2]andl3]l, the uncertainty in the X-ray bore-sight correction, 
and the uncertainty in the astrometry of the /?-band image. The 
resulting Icr uncertainties on the X-ray positions range between 
0".' 1 for the brightest X-ray source to 0'.'7 for the fainter sources. 
In order to identify the optical counterparts to the X-ray sources, 
we treat every star inside the 99% confidence error circle as a 



potential counterpart. These stars are indicated in the CMDs in 
Fig.|4]and Fig.|5]and finding charts are provided in Fig.|6] 

4. Source classification 



As in our previous papers ( Pooley et al. 2002 Bassa et al. 2004 
Ko ng et aL| |2006| l, we first use our astrometry to identify stars 
within the, in this case, 99% confidence circles as possible op- 
tical counterparts to the X-ray sources. We then check whether 
these stars have unusual photometric properties as expected for 
candidate counterparts. Finally we combine the optical photom- 
etry of the candidate counterparts, the X-ray properties of the 
X-ray sources and the combined fx/fopt to classify the proba- 
ble counterparts. These steps are illustrated with Figures l6] l4j l3] 
and It] respectively. For the latter figure, we estimate the visual 
magnitude as V - (B + R)/2 where no V-band magnitudes are 
available. The photometry and the offset between the Chandra 
and optical positions are given for selected optical counterparts 
in Table|5] 

4.1. M55 

Starting with M55, we note that the 99% error circles in Fig. [6] 
include relatively bright stars for CXI, CX7, CX8, CXIO, and 
CX16. The other circles contain no significant stars, or stars that 
we consider too far off-center to be probable counterparts. 

CXI, as remarked above, is the dwarf nova discovered and 



studied by Kaluzny et al. ( 2005 1, whose optical position we find 
to be compatible with the Chandra position, and which we used 
to correct the Chandra coordinates. In Fig. |4] we see that the 
optical counterpart of CXI has strong Ha emission as expected, 
but is less blue than expected for a dwarf nova, being at the blue 
edge of the main sequence. Its X-ray to optical flux ratio is as 
expected for a cataclysmic variable, well above the line log Lx = 
34.0 - O.AMv that roughly separates cataclysmic variables from 
the magnetically active binaries below it (Fig.|7]i. 
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Fig. 5. Optical color-magnitude diagrams for M55 and 
NGC 6366 constructed from the [Anderson et al. (i2008) pho- 
tometry. Candidate optical counterparts to the Chandra X-ray 
sources are indicated and numbered. 



The error circle of CX2 shows a very faint blue object in 
the ACS/WFC observations. The source could be a cataclysmic 
variable, though it is below the detection limit in the FORS2 
observations, so we do not know if it exhibits an excess of Ha 
emission. Compared to classified cataclysmic variables in glob- 
ular clusters, the optical counterpart is significantly fainter The 
probable optical counterpart of CX7 is at a position slightly be- 
low and towards the red of the subgiant branch in the colour 
magnitude diagram F ig. [4} and thus a sub-subgiant in the termi- 
nology introduced by |Belloni et al. ( 1998 1 in their study of M67. 
It has no significant Ha emission. Its fx/fopt indicates a mag- 
netically active binary in agreement with the classification as a 
sub-subgiant. The possible optical counterpart of CX8 from its 
position in the colour magnitude diagram appears to be an ordi- 
nary subgiant; as this star is less well centered in the astrometric 
confidence circle, it is quite possible that the real counterpart of 
CX8 is a fainter object. CXIO has two possible counterparts, the 
brighter one is on the binary sequence in the colour magnitude 
diagram, but puzzlingly has strong Ha absorption. Because the 
R - Ha is not affected by interstellar absorption, this offset from 
the main sequence cannot be explained as due to the star being a 
fore- or background object. This suggests it is an active binary. 
The HST observations show that the single star found in the er- 
ror circle of CXI 6 in the FORS2 observations is in fact a blend 
of three stars. Two of these stars are rather red and too far from 
the main sequence of M55 to be a member. Their red colours 
suggest they are background objects. The third star lies on the 
cluster main-sequence in V - /, and could indicate that the star 
is an active binary. Because the star is blended in the F0RS2 ob- 
servations, it is unclear to which object the Ha absorption can 
be attributed. 

For the sources outside the half-mass radius of M55, CX19 
has two candidate counterparts, one on the subgiant branch, the 
other on the main-sequence. Neither counterpart displays a sig- 



nificant emission or absorption of Ha flux. CX20 on the other 
hand has very strong Ha emission and is very blue in B-R. With 
a X-ray luminosity of Lx = 3.3 x 10^' erg s"' it is very similar 
to CXI and possibly a cataclysmic variable. The counterpart to 
CX21 is on the main-sequence and could be an active binary. 
However, the fx/fopt ratio of this source places it in the region 
that is primarily populated by CVs, the X-ray source could be a 



cataclysmic variable. We will return to this source in Section 5.2 



CX23 shows similarities with CXI 6 in the sense that it is very 
red and has a large excess of Ha emission. It is likely a back- 
ground object. The optical counterpart to CX24 is slightly blue 
but shows no Ha emission or absorption compared to that of the 
cluster stars. It is the brightest X-ray source in our sample, and 
it has an fx/fopt ratio that would be indicative of a cataclysmic 
variable. However, the absence of Ha suggests it belongs to the 
Galactic field. The optical counterpart to CX29 may be yet an- 
other sub-subgiant; it is located below the subgiant branch and 
has a X-ray luminosity comparable to that of CX7 in M55 and 
CX5 in NGC 6366. We do note that the source lies at the edge of 
the error circle. 



4.2. NGC 6366 

Turning now to NGC6366, we consider four possibilities for the 
extended emission of source CXlb; a planetary nebula, a super- 
nova remnant, a group or cluster of galaxies and a nova remnant. 
Planetary nebulae and supernovae tend to have strong Ha emis- 
sion. The apparent absence of such emission in our Ha images 
of CXI (Fig.l2]i argues against these two possibilities. A sec- 
ond argument against a planetary nebula is that the X-ray spec- 
tra of planetary nebula observed with Chandra are rather soft, 
limited to energies below 1.5 keV (e.g. BDh-30°3639, Kastner 
et al."2000', NGC 7027, [Kastner et al.||200T] NGC6543nChu 
et al. 2001), whereas 50% of the photons of CXlb is above 2keV 
(Fig.l2]). A second argument against a supernova remnant is that 
it is an unlikely object to be found in a globular cluster, or as a 
foreground object at a Galactic latitude of 16°. 

The clusters of galaxies studied by Kaastra et al.| ( [2004[ ) and 
[Piffaretti et al.[ (12005 ) have typical X-ray diameters of order 
1 Mpc and X-ray luminosities in excess of 10"^^ erg s"' . Both the 
flux and diameter of CXlb imply a distance of more than a Gpc 
if this source is a typical cluster. The intra-cluster gas in such 
bright clusters strongly dominates the total flux, and so the su- 
perposition of CXI a and CXlb would be accidental if CXlb is 
a galaxy cluster Groups of galaxies have a smaller extent and 
lower X-ray luminosity, against which an individual galaxy can 
stand out (e.g. Nak azawa et al.[2007] l. However, even our deepest 
B and /?-band images show no obvious galaxies near the position 
of CXI, and this excludes a group or cluster of galaxies out to 
several hundred Mpc. 

Finally, a nova outburst may give rise to extended X-ray 
emission, as illustrated by the case of GK Per (Nova Persei 1 90 1 ) 
which was detected as an extended source with Chandra, sur- 
rounding the point source of the cataclysmic variable GKPer 
itself (Balman 2005 1. Both in size and in luminosity CXI is a 
plausible old nova at the distance of NGC 6366, but the spec- 
trum of the extended source CXlb is rather harder than the soft 
extended source surrounding GKPer. The amount of matter in a 
nova remnant is much less than in a planetary nebula, perhaps 
compatible with the absence of extended Ha emission; indeed, 
Slavin et al. ( 1995) l argue that the faint emission around GKPer 
in the Ha filter is due largely to N ii emission. 

For the moment, we consider an old nova to be the more 
likely explanation for the X-ray source CXI in NGC 6366. 
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Novae are known in globular clusters; TSco (1860) in M80, 



V1148Sgr(1943)in NGC 6553 and Nova 1938 in M14 ( Sawyer 



1938 Mayall 1949) Sawyer Hogg & Wehlau 



19641. The nova 



hypothesis explains the position of CXI near the cluster center, 
and the superposition of CXI a on CXlb. It remains to be shown, 
however, that the X-ray spectrum of an old nova can be as hard 
as that of CXlb. 

The two bright possible optical counterparts of CX4 are both 
ordinary main sequence stars according to the colour magnitude 
diagram; the fx/fopt suggests a magnetically active binary, in 
which case the companion does not contribute significantly to 
the optical fluxes. However, CX4 is a hard X-ray source which 
may be hard to reconcile with a magnetically active binary. It 
may be that the source is not related to the globular cluster 
and is assocated with the faint optical counterpart seen in the 
ACS/WFC observations. The optical counterpart of CX5 is yet 
another sub-subgiant as CX7 in M55. For NGC 6366 CX8 co- 
incides with the bright star BD-04°4280. With V = 10.7, it is 
saturated in our FORS2 observations (see Fig.[6| and hence not 
present in Fig.ffland TablelS] At this brightness it is at least 8 mag 
brighter than the cluster turn-off and hence highly unlikely to be 
a member of NGC 6366. Finally CX14 has a candidate counter- 
part that is blue in B-R with respect to the cluster main-sequence 
but has no Ha emission or absorption. The absence of Ha emis- 
sion suggests it belongs to the Galactic field. 



5. Discussion 

5.1. Comparison with earlier observations 

Both M55 and NGC 6366 have been observed with ROSAT, and 
XMM also observed M55. 



5.1.1. M55 



In Fig.lSlwe compare the countrates of sources in M55 observed 
with Chandra (from T able[2 1 with those observed with XMM 
EPIC (from Table 2 of |Webb et aT|j2006). AH XMM sources 



within the field-of-view of Chandra have been detected; the 
XMM sources are the sources with the highest countrates in 
Chandra, plus CX29 (XMM45). Comparison of the Chandra and 
XMM source positions shows that the error estimates by Webb 
et al. for the XMM source positions are overestimated by about 
a factor 2. As expected, the XMM countrate for most sources 
is ~ 4 times higher than the Chandra countrate: for example, a 
2 keV bremsstrahlung spectrum absorbed by the column towards 
M55, with 0.001 cts s"' in Chandra ACIS-S will produce about 
0.0026 cts s"' on the PN and 0.00075 cts s"' on each of the two 
MOS detectors, for a total of 0.0041 cts s"' in XAflMJ^ 

Remarkably, all three sources with the lowest XMM count 
rates were brighter during the Chandra observation. One of 
these is the cataclysmic variable CVl, our source CXI. Webb 
|et al.| ( [2006l l state that this source was three times fainter, m the 
0.5-10 keV band, during the XMM observation than during ear- 
lier ROSAT observations. Since the ROSAT measurements are 
limited to the 0.2-2.5 keV band, we prefer to make the com- 
parison of the luminosities in this energy range. The hardness 
ratio of CXI as measured with Chandra (see Fig.|3]l is that 
of a bremsstrahlung spectrum with temperature kT\, ^ 2keV. 



at 

we 



3 We estimate this with the PIMMS tool 
|http : //heasarc . gsf c . nasa . gov/Tools/w3piimn s.htinl , 
have multiphed the countrates within a 15" extraction radius as given 
by PIMMS with 1 .47 to estimate the total countrates 
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Fig. 7. The 0.5-2.5 keV X-ray luminosity Lx against the abso- 
lute V-band magnitudes. The V-band magnitudes are estimated 
by y = {B + R)/2 to allow a comparison to X-ray sources 
with optical counterparts detected in other globular clusters. The 
X-ray sources are classified into the different classes as indi- 
cated. Besides the Chandra X-ray sources discovered in M55 
and NGC 6366 (which are numbered) we also have included 
classified X-ray sources from 47Tuc, NGC 6752, NGC 288, 
NGC 6397, wCen, M22 and M4. The diagonal dashed line of 
constant X-ray to optical luminosity ratio (logLx - 34.0 - 



0.4Mv) roughly separates CVs from ABs (after Bassa et al. 
^004|). The sohd diagonal line (log Lx = 32.3 - 0.27My) is af- 
ter |Verbunt et al.| ( f2007 ). The two dashed-dotted lines are com- 
puted from 1 1 .2 Gyr isochrones at z = 0.001 and z - 0.019 from 
Girardi et"aL](|2000|), assuming Lx - O.OOlLboi- 



Assuming this spectrum, the distance and absorption of M55 as 
above, and the countrates given in Table 3 of Verbunt| (2001 r, 
cf. Johns ton etaL] ( fT996 ) and in Table 2 of |Webb et aL]p006 i, 
we obtain 0.5-2.5 keV luminosities, in units of 10^' ergs"\ 
of 6.0 ± 1.5 for both the ROSAT PSPC observation (March 
1991) and the ROSAT HRI observation (September 1997), and 
of 0.44 + 0.09 for the XMM observation (October 2001); our 
Chandra observation (May 2004, Table 1) gives 4.1 ± 0.4. Thus, 
the X-ray flux was an order of magnitude lower during the XMM 
observation than during the ROSAT and Chandra observations. 
Webb et al. ( 2006) assume that the X-ray flux is higher dur- 



ing an optical/ultraviolet dwarf nova outburst, and conclude that 
CVl was in outburst during both ROSAT observations. However, 
most dwarf novae have lower X-ray fluxes during outburst than 
in quiescence ( [Wheatley et al.|1996[ Verbunt et al.|1999| and the 



Verbunt et a 



1997 1. We conclude that 



discussion of Fig. 2 in 
CVl was in outburst during the XMM observations, and in qui- 
escence during the ROSAT and Chandra observations, in good 
agreement with the observations by ,Kaluzny et al. ( 2005) that 
CVl spends most of the time in quiescence. 

We detect the optical counterpart to CVl at R - 20.7 and 
B- R - 1.34, where it lies only about 0.3 mag blue-wards of the 
cluster main-sequence. With R — Ha = 0.58, it exhibits a large 
excess of Ha emission, reminiscent of cataclysmic variables. At 
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Table 5. Candidate optical counterparts to the Chandra X-ray sources. The positional offsets are given relative to the Chandra 
values in Tablesl2]and[3] The X-ray luminosity is in the 0.5-2.5 keV energy range. The first optical counterpart to CX16 is a blend 
of the three other stars. 
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B - 22.1, our photometry of CVl puts it close to the quiescent 
level of about B 



-2.5 -2 

log counts/s (XMM) 

Fig. 8. Detected countrate in the 0.5-6.0 keV energy band (Xsoft+ 
Xhard, Table[2]) of the Chandra ACIS-S observation as a function 
of the XMM EPIC (PN and both M PS detectors) countrate in the 
0.5-10 keV energy band (Table 2 of |Webb et al .T2006 ). The soHd 
and dashed lines show equal countrates and 4 times higher/lower 
countrates, respectively. For most sources, the XMM countrate is 
about 4 times higher than the Chandra countrate, as expected for 
a constant source. 



22.5 seen in the light-curve by Kaluzny et al. 
( |2005) l. Furthermore, the proximity in color of the source to the 
cluster main-sequence was also seen by Kaluzny et al. (2005i 
during quiescence, where the source actually coincides with the 
main-sequence in B - V and V-I (see Fig. 4 of Kaluzny et 
al.). Hence, we conclude that CXl/CVl was also in quiescence 
during our FORS2 observations. 

Fig.|8] shows that CX20/XMM14 had Chandra and XMM 
countrates as expected for a constant luminosity between the two 
observations. The source was not detected with ROSAT; in re- 
analysing the ROSAT PSPC observation, using the method out- 
lined inlVerbunt (2 001| l, we obtain a 2cr upper limit to the coun- 
trate of 0.0006 ctss" (channels 50-240), which corresponds 
roughly to a flux at 0.5-2.5 keV of 8.9 x 10"'^ ergcm"^ s"'. 
Analysis of the HRI data gives an upper limit of 10 source 
counts, or 8 X 10"'^ergcm"^s"'. Within the error, these limits 
are compatible with the flux found with Chandra (Tablel2]i and 
with XMM. |Webb et aL] ( |2006p argue that this source has been 
variable by at least a factor 22 between the ROSAT and XMM 
observations. We find no evidence for such variation. 



[Kaluzny et al.| ( |2005| l discuss an interesting blue variable 
star, which they name M55-B1. From its variability between 
observing seasons, in the absence of short-term (~ hr) variabil- 
ity, they argue that it is a background quasar rather than a cat- 
aclysmic variable. This source is not detected on the Chandra 
ACIS-S3 chip, but on an adjacent chip. It is clearly detected in 
X-rays, with an unabsorbed flux in the 0.5-8.0keV band of 5.1 x 
10"'^ ergcm"^s"' for absorption fixed at that of M55. A free 
power law fit gives a higher absorption (A^h = 7.4 x 10^' cm"^) 
and higher unabsorbed flux of 8.8 x 10"'^ erg s"' . 
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5.1.2. NGC 6366 

Only the brightest Chandra source was detected with ROSAT, at 
a 0.5-2.5 keV luminosity f Verbuntl |2001) l compatible with that 
during the Chandra observation. ROSAT observed CXI a and 
CXlb as one source. The observed spectrum of the point source 
is relatively soft (see Fig.|2] and Sect. 2.3 i, which with the high 
absorption towards NGC 6366 implies a very soft intrinsic spec- 
trum. Whereas many cataclysmic variables have hard spectra, a 
very soft spectrum is still compatible with identification as an 
old nova. 



5.2. Ratio of optical to X-ray luminosity 

Our classification of the X-ray sources for which we have an 
optical identification is based in part on the X-ray to optical lu- 
minosity ratio. |Bassa et aL] ( |2004| showed that a line of constant 
Chandra X-ray to optical luminosity roughly separates the cata- 
clysmic variables from the active binaries in 47 Tuc, NGC 6752 
and M4 (data from JGrmdlay et al.|200T] [Edmonds et al.|2003 



Pooley et al. 2002 Bassa et al. 2004 1. This separatrix corre- 
sponds to the dashed line in Fig. 171 and is given by 



log Lx(erg s"') [0.5-2.5 keV] = 34.0 - 0.4Mv 



(1) 



Remarkably, and perhaps disconcertingly, this line is rather 
higher than the upper bound to the X-ray luminosities of stars 
and RS CVn binaries near the Sun. An upper bound is given in 
IVerbunt et aL] ( |2007l i as 



logLx(ergs~i) [0.5-2.5 keV] = 32.3 - 0.27Mv 



(2) 



This line is an absolute upper bound in that it lies above all 
ROSAT measurements for nearby stars (selected from Hiinsch 
eraL][T999l l and for RS CVn systems ( [Dempsey et al.| |1993| l, 
where we multiply the 0. 1-2.4 keV luminosities given by 
Hiinsch et al. and Dempsey et al. with 0.4 to obtain the 0.5- 
2.5 keV luminosities. It is shown as a solid line in Fig. [7] It has 
been remarked already on the basis of Einstein 0.2^ ke V mea- 
surements that the X-ray luminosity of main-sequence stars in- 
creases with increasing stellar rotation velocity, until it saturates 
at 



Lx - O.OOlLboi 



(3) 



where Lboi is the bolometric luminosity of the star (e.g. Vilhu & 
Walter|1987| esp. their Fig. 6). 

We show in Fig.[7] th e line of O.OO lLboi for the 11.2Gyr 
isochrones calculated by [Girardi et al.| (|2000) for metal -poor 
stars (z = 0.001), appropriate for globular clusters, and for stars 
of solar abundance (z - 0.019), appropriate for stars in the solar 
neighbourhood. These lines differ little, and give a more conser- 
vative upper bound to the maximum 0.5-2.5 keV luminosity of 
stars near the Sun. 

The region between the separatrix Eq.fTlon the one hand and 
the upper bounds given by Eq.l2]or Eq.l3]contains a dozen opti- 
cal counterparts of X-ray sources in globular clusters that have 
been classified as magnetically active binaries. If this classifica- 
tion is correct, we must accept that such binaries can have rather 
higher X-ray luminosities in globular clusters than in the solar 
neighbourhood. Since the bolometric luminosity as a function of 
My is similar for low and solar metallicities, this would imply 
that the saturation of X-rays occurs at higher Lx/Lboi in globular 
clusters. Alternatively, we may wish to reconsider the classifica- 
tion of these sources. If we want to classify them as cataclysmic 
variables, we must accept that globular clusters host cataclysmic 



Table 6. Values for central density po, core radius r^, distance d, 
and absolute visual magnitude My, taken from Harris (1996, ver- 
sion Feb 2003), and for M4 from |Richer eraL] ( [l997 ). ColHsion 
number F and mass within the half-mass radius Mb are scaled on 
the values for M4. 



cluster logpo ''c d 

iU pc-') (") (kpc) 



My 



Mh 



M4 


4.01 


49.8 


1.73 


-6.9 


= 1.00 


= 1.00 


NGC 6397 


5.68 


3.0 


2.3 


-6.6 


2.05 


0.78 


47 Tuc 


4.81 


24.0 


4.5 


-9.4 


24.91 


10.19 


NGC 288 


1.80 


85.2 


8.8 


-6.7 


0.04 


0.86 


M55 


2.15 


169.8 


5.3 


-7.6 


0.18 


1.82 


NGC 6366 


2.42 


109.8 


3.6 


-5.8 


0.09 


0.35 



variables that are not blue with respect to the main sequence, 
even in the B-band filter. These systems clearly warrant further 
study. We note that CVl in M55 does lie on the main-sequence 
in B - y in quiescence (Kaluzny et al.|20()5 see Sect. 5.1.1 1. 



5.3. Origin of the X-ray sources 

To compare the numbers of X-ray sources in globular clus- 
ters, Pooley et al.^ (2003) chose a luminosity limit of Lx > 
4 X 10"" erg s~' [0.5-6.0 keV], which is detectable in almost all 
the clusters they investigated. We have detected 16 sources above 
this limit within the half-mass radius of M55, with an expected 
number of 8-9 unrelated sources; and 5 within the half-mass 
radius of NGC 6366, with an expected number of 4 unrelated 
sources. As akeady argued above, from these numbers alone 
we cannot exclude that all detected sources are unrelated. The 
situation changes if we take into account our optical identifica- 
tions (Section 4), which indicate that at least three sources (CXI, 
CX7, CXIO) are highly probable members of M55, and that two 
sources (CX4, CX5) are highly probable members of NGC 6366. 
We have also argued that the point source plus extended source 
is an old nova in NGC 6366. 

In Table|6]we compare the collision numbers F = po'^'"c^ 
and the masses within the half-mass radii M^ (as measured with 
the absolute magnitude assuming a fixed mass-to-light ratio) for 
several well-observed clusters; where we scale both F and M^ on 
the globular cluster M4. 

: 2 sources above the luminos- 



Pooley et al. (2003 1 report 41 



ity limit in the core of 47 Tuc. If the number of X-ray sources 
would scale with the collision number, we would predict 0.3 
sources in M55, and 0. 15 in NGC 6366. The number of detected 
very probable member X-ray sources is significantly higher in 
both clusters. If the number of X-ray sources would scale with 
the mass, we would predict about 7 sources in M55, and 3 to 4 
in NGC 6366. The number of detected member X-ray sources is 
remarkably close to this. 

To study the origin of X-ray sources in globular clusters in 
further detail, [Pooley & Hut| ( |2006[ ) investigate the hypothesis 
that the number of sources scales with the collision number and 



N = ar + bM 



(4) 



and show that the dependence on mass is significantly present 
in the observations. To overcome the small-number statistics of 
clusters with small source numbers, they combine the numbers 
of such clusters before fitting them to Eq.H] 

In a recent paper, Verbunt et al. (2007 ) suggest a more gen- 
eral way of fitting the observed numbers to Eq.|4] Briefly, they 
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determine for each cluster the most probable combination of 
numbers N^ of member and A^b of unrelated sources that gives 
the observed number of sources within the half-mass radius A^h, 
based on the expected number of cluster sources fi^ according 
to the model of Eq.H] and on the expected number of unrelated 
sources /^b- In doing so they assume that the probability of find- 
ing A^ sources when fi are expected is given by the Poisson func- 
tion 



PiN,^.) = ^.- 



(5) 



The fitting procedure consists of varying a and b to maximize 

P^Yl^P(N,,^i,)P{Nb,fib)]j (6) 

7 

where j indexes the clusters. 

The values for A^b and //(. that we use are taken from |Pooley 
et al. ( 2003| l. We ignore the clusters for which only a lower 
limit to A^h is known. We add three clusters: NGC 288, M55 and 
NGC 6366; the values we use are 11/8.0, 16/8.3 and 5/4.0 for 
Nhluc for these three clusters, respectively. The best solution is 
found for a - 1.3rM4 and b - 0.6Mm4. The contours of the ac- 
ceptable solutions are given in Fig. 9\ and it is seen that /? = is 
acceptable at about the 2cr level. Thus, from the observed num- 
bers alone, we have no firm evidence for a dependence on mass. 

We repeated the fitting procedure, but excluding solutions 
for which A^c is less than the minimum number of highly proba- 
ble cluster members. For these fits we assume that 2 sources are 
certain members of NGC 288, 3 of M55 and 2 of NGC 6366. As 
expected, these constraints force the solutions to a stronger de- 
pendence on mass: the best solution is found for a - 1.2rM4 and 
b = 1.1Mm4. As seen in Fig.l9j the dependence on mass is now 
significant. 

The best values of a and b are relatively stable against small 
variations in the numbers jj^ and A^h for individual clusters, as 
comparison with the results of Verbunt et al. ( 2007 | l confirm. The 
probability of the best solution does vary with small variations 
in jjc and A^h for individual clusters. For example, increasing the 
minimum of cluster member X-ray sources in NGC 6366 from 2 
to 3 does not change the best values for a and b, but decreases 
the probability of the best solution for NGC 6366, for which the 
best model predicts jj^ - 0.4, and thus also slightly decreases the 
overall probability. 

The cluster which is responsible for the lowest probability is 
NGC 6397, which houses 12 sources, of which our best model 
assigns A^c = 11 to be cluster members, even though it expects 
only jUc = 3.6. The remarkably large number of sources in this 
cluster has been remarked upon before, and various explanations 
have been offered. In our model, however, whereas still the worst 
fit cluster, it is not altogether unacceptable. The main difference 
between this result and the strong deviation for NGC 6397 as 
found by e.g. |Pooley et al.|(; 2003) appears to be due to a differ- 
ent definition of F. Pooley et al., {2003) integrate a King model 



to estimate F, whereas we compute it as F = po 'c ■ Our prob 
abilities do not take into account the uncertainties in F and Mb, 
nor the observed number A^b which through its dependence of 
the flux limit is affected by the assumed distance to the cluster 
Clearly these uncertainties warrant further investigation. 

An important result illustrated by Fig.l9]is that the masses 
of the clusters may not be taken as proxies for their collision 
numbers. Even though there is some correlation between mass 
and collision number, the spread in collision number at given 
mass, and in mass at given collision number, is very large. Indeed 




Fig. 9. Probability contours around the best solutions for a fit of 
the observed number of sources within the half-mass radii of 12 
globular clusters to Eq.fflwhere we allow zero cluster members 
in some clusters (dashed contours), and where we impose lower 
limits to the number of cluster members of 2 for NGC 288, 3 
for M55 and 2 for NGC 6366 (solid contours). We assume that 
-2{P - Pbf), with the probability P given by Eq.l6]and Pbf indi- 
cating the best fit probability, follows ax^ distribution, and show 
Ax^ = 1,4 and 9. The imposition of minimum numbers of clus- 
ter member X-ray sources strongly increases the significance of 
the mass dependence in Eq.|4] 



our fits show that the number of X-ray sources within a cluster 
is very badly predicted by the mass alone, as illustrated by the 
highly significant offset of the best solution from the a -0 line. 
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Fig. 6. Finding charts for all Chandra X-ray sources in M55 and NGC 6366 which coincide with the F0RS2 or ACS/WFC field-of- 
view (see Fig.fT]). The first 26 charts, CXI to CX29 are for sources in M55 and are 8" x 8" in size. Sources CXI to CX17 are located 
in the cluster half-mass radius, the other sources are located outside of it. The next 14 charts are for all X-ray sources in NGC 6366 
and these are 10" x 10" in size. The first 5 are for sources inside the half-mass radius, the other 9 outside the half-mass radius. 
The ACS/WFC images have a 5 times better pixel scale and show the drizzled F606W observations. The finding charts for the other 
sources were constructed from the R-hand FORS2 observations. The uncertainties on the X-ray source position is indicated with a 
99% confidence error circle. All images have North to the top and East to the left. 



